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I.  INTRODUCTION 


The  reaction  zone  in  cw  chemical  lasers  [Fig.  1(a)]  is  generally 
maintained  at  pressures  of  the  order  of  1  to  10  Torr  to  permit  fast  mixing  of 
the  reactants.  At  these  pressures,  the  spectral  lineshape  is  Inhomogeneous ly 
broadened,  i.e.,  the  radiation  field  interacts  with  only  a  portion  of  the 
excited  molecules.  As  a  result,  hole  burning*  may  occur  as  the  degree  of 
optical  saturation  is  Increased.  Hole  burning  affects  laser  output  power  and 
the  index  of  refraction  of  the  lasing  medium. 

A  comprehensive  theory  for  inhomogeneous  broadening  effects  in  a  steady- 
state  laser  oscillator  has  been  developed  by  Lamb.  Lamb's  theory  was  gener¬ 
alized  in  Ref.  3  to  account  for  effects  of  cross  relaxation  and  streamwlse 
flow  variations  in  cw  chemical  lasers.  Results  were  obtained  for  the  case  of 
a  Fabry-Perot  (F-P)  resonator  with  a  single  longitudinal  mode.  The  single 
longitudinal  mode  case  corresponds  to  a  mirror  separation  distance  L 
[(Fig.  1(a)]  up  to  approximately  1  m.  For  high-power  lasers,  the  mirror 
separation  can  be  of  the  order  of  10  m,  and  many  longitudinal  modes  are 

excited.  The  latter  case  was  treated  in  Ref.  4  in  the  limit  Av  «  Av,  , 

c  n 

Av^  «  Av^  where  Av^  and  Av^  are  characteristic  homogeneous  and  Doppler 
widths,  respectively,  and  Avc  is  the  longitudinal  mode  spacing  (Fig.  2).  In 
Refs.  3  and  4,  a  simplified  two-level  model  is  used,  and  closed  form  solutions 
are  obtained.  A  numerical  code  developed  by  Bullock  and  Lipkis^  treats  the 
full  system  of  equations  that  describe  inhomogeneous  broadening  effects  in  cw 
chemical  lasers.  Results  from  Refs.  4  and  5  are  in  agreement  in  those  regimes 
where  both  are  applicable.^ 

Transverse  flow  expansion  (i.e.,  source  flow)  can  be  used  in  cw  chemical 
lasers^  to  reduce  the  temperature  increase  in  the  lasing  region.  The  result¬ 
ing  mean  motion  in  the  optical  path  direction  (e.g..  Fig.  3)  modifies  the 
spectral  lineshape  of  the  lasing  medium.  The  effect  of  source  flow  on  a 
Doppler  broadened  lineshape  is  deduced  in  Ref.  6  for  the  case  of  small  lateral 
motion  and  in  Ref.  7  for  the  case  of  arbitrary  lateral  motion. 


Figure  2(b).  Characteristic  frequencies  and  frequency-dependent 

variables  for  case  Av  «  Av  ,  Av  «  Av  ;  lasing 

tl  a  C  n 


n 

i 
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Figure  3.  Source  flow 
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In  this  report,  the  multiple  longitudinal  mode  F-P  resonator  theory  of 
Ref.  4  ia  modified  to  account  for  source  flow  effects.  The  primary  objective 
is  to  evaluate  the  power  decrement  associated  with  source  flow  in  cw  chemical 
lasers.  Numerical  results  are  presented  for  cw  chemical  lasers  with  laminar 
diffusion. 
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II.  THEORY 


The  model  of  Refs*  3  and  4  Is  generalized  herein  to  Include  non- 
Maxwe Ilian  Doppler  lineshapes.  The  performance  of  a  multiple  longitudinal 
mode  cw  chemical  laser.  Including  source  flow  effects,  is  then  deduced.  The 
development  follows  that  in  Ref.  4.  Unless  otherwise  noted,  the  notation  is 
the  same  as  used  in  Ref.  4.  Symbols  are  defined  in  Appendix  A. 


A. 

REGIME 

It  is  assumed  that 

Av.  <<  Av , 
n  a 

(la) 

Av  <<  Av. 
c  h 

(lb) 

where  Av^  and  Av^  are  characteristic  homogeneous  and  Doppler  widths,  respec¬ 
tively,  and  Avc  is  the  longitudinal  mode  spacing,  as  illustrated  in  Figs. 

2(a)  and  (2b).  Expressions  for  these  quantities  are  presented  in  Ref.  3.  The 
inequality  in  Eq.  (la)  is  generally  satisfied  in  cw  chemical  lasers  and 
simplifies  the  relation  between  small  signal  gain  and  particle  number  density 

O 

distribution.  Equation  (lb)  is  used  to  simplify  the  effect  of  F-P  resonator 
modes  on  lineshape.^  The  present  results  are  believed  valid^  for 
Avc/Av^  <  0(1).  The  latter  inequality  is  satisfied  for  F-P  mirror  separation 
of  the  order  of  L  >  0(10)m. 

B.  MODEL 

A  cw  HF  chemical  laser  is  illustrated  in  Fig.  1(a).  The  flame  sheet 
model  of  Refs.  3  and  4  is  illustrated  in  Fig.  1(b).  The  reactants  are  assumed 
to  be  premixed  but  do  not  react  until  a  flame  sheet  y^(x)  is  reached.  The 
form  of  the  flame  sheet  is  specified,  a  priori,  from  diffusion  theory.  The 
streamwise  station  where  the  flame  sheet  reaches  the  channel  centerline, 
denoted  Xq,  characterizes  the  diffusion  rate. 


The  reactants  (e.g. ,  H2  +  F)  form  upper  vibrational  level  species  at  the 
flame  sheet.  Characteristic  rates  are  normalized  by  the  collisional  deactiva¬ 
tion  rate  kc<j.  Streamwise  distance  is  expressed  in  the  form  £  ■  k£d  x/u, 
which  is  the  ratio  of  a  convection  time  to  a  collisional  deactivation  time  and 
is  an  order  1  quantity.  The  quantity  £^  ■  kfid  x^/u  Is  the  ratio  of  the  diffu¬ 
sion  time  to  the  collisional  deactivation  time.  It  is  assumed  that  £^  >  £fi, 
where  £g  denotes  the  station  at  which  lasing  is  ended. 

The  difference  between  the  net  population,  per  unit  volume,  in  the  upper 
and  lower  lasing  levels  is  expressed  in  normalized  form  as 

AN  =  N2  -  N1  -  (n2  -  a1)  yf/(nrw)  (2a) 

where  nr  is  a  characteristic  reactant  number  density  upstream  of  the  flame 
sheet,  e.g.,  F  in  Fig.  1(b).  Particles  (per  unit  volume  per  unit  frequency) 
resonant  with  laser  frequencies  in  the  range  v  to  v  +  dv  are  denoted  n(v)  and 
are  normalized  in  the  form 

ANy  =  N2(v)  -  NL(v)  »  ln2(v)  -  n^v)]  yf/(nr‘p()  w)  (2b) 


where 

pQ  -  [4  (In  2)/it]1/2/Avd  (2c) 

The  variation  of  AN  with  v  is  illustrated  by  Figs.  2(a)  and  2(b)  for  non- 

v  4 

lasing  and  lasing  cases,  respectively. 

The  quantities  AN  and  ANy  are  related  by 


AN 


f  ANV  ^ 


(3) 
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The  normalized  gain  per  unit  length  G(v)  equals 


G(v) 


g(v)  yf 


a  n  w  p  Av. 
o  r  ro  h 


ANv> 


(4a) 


where  £(v-v' )  Is  the  Lorentzian  (homogeneous  llneshape) 


5?(v-v') 


g(v,v<  ) 

o 

o 


(4b) 


The  quantity  g(v)  In  Eq.  (4a)  Is  the  gain  per  unit  length  In  the  reactive  flow 
region  0  <  y  <  y^  [Fig.  1(b)].  The  average  value  of  gain  per  unit  length. 
Including  both  reactive  and  nonreactive  flow  regions,  is 

g(v)av  -  g(v)  yf/w  (4c) 

Thus,  the  quantity  G(v)  is  seen  to  be  a  normalization  of  the  average  gain  per 
unit  length. 

At  this  point,  we  depart  from  the  development  in  Ref.  4  and  assume  a 
source  flow  with  a  linear  variation  of  transverse  velocities  in  the  range 
-  V@  <  V  <  where  Vg  is  the  transverse  velocity  at  the  edge  of  the  source 
flow  (Fig.  3).  In  the  absence  of  radiation,  the  particles  have  a  distribu¬ 
tion^ 

ANv 

“Jg-  "  F(X,X@)  (5a) 


where 


F(X,Xg) 


[erf  (X  +  X  )  -  erf  (X  -  X  )] 
e  e 


(5b) 
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(5c) 


-  e  (X  -  0) 

e 

-  (»l/2/2X  )  erf  X  (X  -  0) 

e  e 

X  -  2  (In  2)1/2  (v  -  v  )/Av  , 

o  a 


(5d) 

(5e) 

(5f ) 


The  function  F(X,Xe)  is  normalized  such  that 

(2/tt1/2)  J  F(X,X0)  dX  -  1 


(5g) 


The  quantity  Xg  is  the  ratio  of  maximum  transverse  velocity  Vg  to  the  most 
probable  particle  thermal  speed  a  and  is  a  measure  of  the  importance  of  the 
transverse  motion  relative  to  thermal  motion.  Note  that  F(X,Xfi)  is  Maxwellian 
when  Xe  ■  0.  In  the  Doppler  approximation  (Av^  «  Av^)  ,  the  zero  power  gain 
equals 

G(v)  -  (*/ 2)  ANy  -  (ir/2)  AN  F(X,Xg)  (6) 


which  is  illustrated  in  Fig.  4(a).  Line  center  gain  is  decreased  and  fre¬ 
quency  width  is  increased  as  X.  increases.  If  we  let  v.  denote  the  center 
frequency  for  each  longitudinal  mode  in  an  F-P  resonator  [Fig.  2(b)],  the  gain 
equal  loss  condition,  for  each  mode,  can  be  expressed^ 


G(Vj)  5  Gc 


(-1)  (In  R  )/(w  n  ) 
m  sc 

°o  nr  Avh  P0 


(7a) 
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Figure  4(a).  Variation  of  number  density  with  frequency  at 
two  streamwise  stations;  upstream  of  lasing 
region,  C  <  C, 


where  Is  the  reflectivity  of  one  mirror  and  nflCw  is  the  net  width  of  the 
gain  region  [Fig.  1(a)].  Equation  (7a)  follows  from  the  average  value  of 
gain,  under  threshold  conditions  in  an  F-P  resonator,  i.e., 

(8c)av  "  t-M**  Ra)/(wnsc)  (7b) 

For  convenience,  the  mirror  separation  L  is  assumed  to  permit  a  longitudinal 
mode  at  vq.  Lasing  is  initiated  at  the  streamwise  station  where  G(vQ)  ■  Gc. 
If  we  let  AN^  denote  the  value  of  AN  and  ANy  denote  the  line  center  (X  *  0) 
value  of  ANy  at  the  station  where  lasing  is  initiated,  these  quantities  are 
related  by 


(2/if)  Gc  -  ANy  -  ANt  F(0,Xe)  (8) 

Equation  (8)  defines  upstream  boundary  conditions  for  the  lasing  region. 

In  view  of  Eq.  (lb),  the  spectral  lineshape  under  laser  conditions  can  be 
assumed  to  have  the  form  [Fig.  4(b)] 

-  1  |X|  <  Xf  (9a) 


F(X,X  ) 

1x1  >x£ 


(9b) 


where  Xf  is  the  value  of  X  corresponding  to  the  largest  lasing  frequency  and 
I X |  <  Xf  denotes  the  lasing  region.  In  Eq.  (9a),  small  departures  from  1  in 
the  lasing  region  are  neglected  [e.g. ,  Figs.  2(b)  and  4(b)]. 
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If  we  let 


H(Xf ,Xfi) 


(10a) 


substitution  of  Eq.  (9)  into  Eq.  (3)  yields 


H(Xf ,X  ) 


1 

F(Xf,Xe) 


F(X,Xe>  dX 


(10b) 


2  Y 

T7T  xf 


+ 


1  * 

1 

■(VXe)2  "(Xf+Xe)2' 

2X  F(X.,X  )  i 

T7I 

6 

e  f*  e  | 

ir 

+  (X,+x  )  erfc  (X,+X  )  -  (X,-X  )  erfc  (X  -X  ) 
fe  fe  fe  fe 

X2 

-  ~fl  xf  +  e  f  erfc  Xf  (Xg-  0) 


-  l/F(0,Xg)  (Xf  -  0) 


(10c) 


(lOd) 


(lOe) 


The  quantities  Xf  and  AN  are  functions  of  streamwise  distance,  which  are 
evaluated  in  the  course  of  the  solution. 

C.  LASER  PERFORMANCE 

The  variation  of  AN  and  AN  with  normalized  streamwise  distance  c  was 

V  -x2 

deduced  in  Ref.  4  for  the  case  of  a  Maxwellian  zero  power  lineshape  e  . 

_y2 

These  equations  are  applicable  in  the  present  study  if  e  is  replaced  by 
F(X,Xe).  Equations  (12a)  and  (12b)  of  Ref.  4  then  become 
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dANv 


F(X,Xe)  [(dUp/dc)  -  Up  +  RAN] 


-  (1+R)  AN  -  ir  Av.  AN  I(v) 

v  n  v 


(lla) 


dAN.fi 
d?  d? 


N„ 


where 


(lib) 


dP  _  2 

dc  =  ~U2 

n 


Avh 


G  I  dX 
c 


(11c) 


Effects  of  chemical  pumping,  collisional  deactivation,  cross  relaxation 
(velocity  exchanging  collisions)  and  stimulated  emission  and  absorption 
as  discussed  in  Ref.  4  are  included  in  Eqs.  (lla)  through  (11c).  Here 
Up  -  Np  +  N2,  R  ■  kcr/kcd  ■  cross  relaxation  rate/collisional  deactivation 
rate,  P  ■  normalized  output  power  extracted  up  to  station  ?,  and  I(v)  repre¬ 
sents  a  continuous  intensity  distribution,  which  replaces  the  discrete  inten¬ 
sity  distribution  Ij  [i.e.,  I(v)  ■  I^/Av^  in  accord  with  the  approximation 

Av  «  Av. .  The  cross  relaxation  parameter  R,  which  can  be  viewed  as  the 
c  h 

number  of  collisions  required  to  deactivate  an  excited  particle,  is  of  the 
order  10  to  100  for  an  HF  laser. ^ 

The  variation  of  X^  with  C  can  be  found  by  consideration  of  Eq.  (lla)  in 
the  region  |X|  >  Xp,  1*0.  Substitution  of  Eq.  (5b)  into  Eq.  (lla),  with 
1*0,  yields 


(-1)  dF(Xf,Xe)/dX^ 
tF(Xf,Xe)]2 


-  (1+R? 

F*Xp^ 


(12a) 
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where  H  is  as  defined  in  Eq.  (10b)  and,  from  Eq.  (5b) 


dF(Xf,Xe) 


(-n 

4  X  X. 
e  f 


'(XfV 


-(Xf+xe)i 


The  initial  condition  for  Eq.  (12a)  is 

Xf  -  0  at  5  -  (12b) 


where  Is  Che  station  at  which  lasing  is  initiated.  The  use  of  Xf  as  the 
dependent  variable  in  Eq.  (12a)  avoids  singular  behavior  near  Equations 

(12a)  and  (12b)  can  be  integrated  to  provide  X^  as  a  function  of  €  if 
Hp/AN^  is  specified.  The  magnitude  of  Xf  increases,  reaches  a  maximum,  and 
then  decreases  as  ?  is  increased.  The  corresponding  variation  of  AN /AN 
with  c  is  then  found  from  Eq.  (10). 


The  lasing  Intensity  I  is  evaluated  by  consideration  of  Eq.  (11a)  in  the 
region  |X|  <  X^,  dANy/dc  ■  0.  Substitution  of  the  latter  into  Eq.  (11a) 
yields 


IT  Avh  I  +  1  +  R 

F(X,X  ) 
e 


+  RH(X.,X  ) 
r  e 


(13a) 


The  right  hand  side  of  Eq.  (13a)  is  Independent  of  X.  Hence,  at  each  stream- 
wise  station,  v  Av^  I  is  equal  to  a  term  proportional  to  F(X,Xe)  minus  the 
quantity  (1  +  R).  At  some  streamwlse  stations,  I  is  negative  in  the  vicinity 

a# 

of  X  -  Xf.  Conditions  under  which  I  is  negative  can  be  deduced  from 
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*  Av,  I 

tl 

F(X,Xe) 


(-1) 

[F(Xf,Xe)l2 


dF(Xf,Xg)  dX2 


F(X,Xe)  -  F(Xf,Xe>' 
F(X,Xe)  F(Xf,Xg) 


(13b) 


2 

which  follows  from  Eqs.  (la)  and  (13a).  The  coefficients  of  dX^/d?  and  1  +  R 

in  Eq.  ( 13b)  are  positive  quantities.  Hence  I  is  positive  for  all  |X|  <  Xf 

2  2  1 
when  dX^/dc  is  nonnegative.  At  downstream  stations  where  dX^/d?  is  negative, 

I  becones  negative  in  the  vicinity  of  X  =  Xp  The  extent  of  the  negative  I 

region  depends  on  the  magnitude  of  R  (i.e.,  the  larger  the  magnitude  of  R  is, 

the  smaller  is  the  extent  of  the  negative  region). 

Negative  values  of  I  correspond  to  power  absorption  by  the  lasing  medium 

2 

[Eq.  (11c)].  The  absorption  is  needed  at  downstream  stations  (dX^/dc  <  0) 
to  maintain  the  constant  gain  boundary  condition  [Eq.  (7a)]  for  longitudinal 
modes  near  X  =*  Xp  Negative  values  of  I  are  unrealistic  and  arise  at  the 
downstream  stations  because  of  the  a  priori  choice  of  llneshape  in  the  region 
I X |  >  Xp  The  present  solution  is  valid  for  those  cases  where  negative  values 
of  I  are  absent  or  negligible.  The  net  power  liberated  up  to  station  c  is, 
from  Eq.  (lib) 


2P 

AN 

v 


H(Xf,Xe) 


+  H(Xf,Xg) 


(14) 


A  similar  expression,  deduced  from  Eqs.  (11c)  and  (13),  indicates  the  self 
consistency  of  the  present  solution.  The  total  output  power  from  the  laser  Pg 
is  found  by  evaluation  of  the  upper  limit  in  Eq.  (14)  at  the  station  where 
dP/d?  -  0  [Eq.  (lib)]. 

The  performance  of  a  cw  chemical  laser,  including  source  flow  effects,  is 
defined  by  Eqs.  (12)  through  (14). 
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D.  LAMINAR  DIFFUSION 

We  now  consider  Che  case  In  which  Che  flame  sheeC  in  Fig.  1(b) 
corresponds  Co  laminar  diffusion.  Arbitrary  values  of  Che  paramecers  R  and  Gc 
are  considered,  and  analytic  solucions  are  Chen  obcalned  for  limicing  cases. 

For  laminar  diffusion,  Che  flame  sheec  has  Che  form 

yf/w  -  (x/xD)^2  (15) 


and  NT  is  approximaCed  by 


3,4 


(16) 


The  race  of  chemical  pumping  is  defined  by  Eq.  (16),  in  which  ic  is  assumed 
chaC  Che  pumping  reacCion  goes  Co  completion  at  Che  flame  sheet. 

Equations  (12)  through  (14)  can  now  be  solved. 

1.  ZERO  POWER  (5*5) 

Integration  of  Eq.  (lib)  in  Che  absence  of  lasing  and  substitution  of 
Eq.  (6)  yields 

Gp/2  AN  -  (,2/it )  ?p/2  G(v)/F(X,Xe)  -  2D(c1/2)  -  <-1/2  (17a) 


where  D  (  )  is  Che  Dawson  integral 


D(X) 


-X* 


/ 


(17b) 
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In  the  absence  of  lasing,  the  net  inversion  AN  Is  Independent  of  Xfi ,  whereas 
the  gain  G(v)  depends  on  Xg.  The  quantity  AN  reaches  a  maximum  at  c  ■  0.3051. 
At  this  station, 

'i/2  ’  °-3528  (18> 

where  subscript  mzp  denotes  maximum  zero  power  value.  Equation  (18)  can  be 
used  to  convert  normalized  variables  to  physical  variables  for  cases  where  the 
maximum  zero  power  value  of  average  line  center  gain  for  the  case  X  ■  Q, 
which  is  denoted  by  has  been  evaluated  numerically  or  experimentally. 

Then,  from  Eqs.  (4a)  and  (18) 


—  1/2 

a  n  p  Av. /zt  m  1.804  g 
o  r  ro  h  D  6mzp 


(19a) 


(« 

c  av  mzp 


1.804  c 


1/2 


(19b) 


where  1.804  ■  (2/ir )/0.3528  and 

®mzp  '  j  -  S -v0^av,  X  “  ol 

i  e  Jmzp 

Equations  (19a)  and  (19b)  apply  for  the  case  of  a  laminar  diffusion  flame  if 
CD  >  0.3051. 

The  zero  power  gain  decreases  to  zero  at  ?  ■  1.1301,  so  lasing  is 
restricted  to  the  region  0  <  z  <  1.1301.  Lasing  is  initiated  at  the  station 
where  the  threshold  gain  is  reached.  Thus, 
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4/2  ANi  -  Gc  -  2D(cJ/2)  -  cj/2 


(20a) 


where 


G 

c 


2 

r 


.1/2 


F(0,Xe) 


(20b) 


For  a  given  value  of  G^,  is  obtained  from  Eq.  (20a)  by  iteration  and  must 

lie  in  the  range  0  <  <  0.3051.  Note,  from  Eqs.  (18)  and  (20),  that 

G  <  0.3528.  The  laser  is  saturated  when  G  /0.3528  «  1.  In  this  region 
c  c 


1/2  '‘w  e  \ 

H  “  GC  1  +  (4/3)  <  +  0 


(20c) 


2.  POWER  ON  (C  >  Ct) 

The  variation  of  Xf  with  ?  in  the  lasing  region  5  >  Cj  is  found  by  sub¬ 
stitution  of  Eq.  (16)  into  Eq.  (12).  The  result  is 


[F(Xf,Xe)]2 
(-1)  dF/dX2 


1  -  2Z 


F(0,Xe) 


+  2ZRH(Xf,Xg)  - 


2Z(1+R) 

wprj 


(21a) 


with  boundary  conditions 


Xf  -  0  at  Z  -  Z  (21b) 

1/2  2 

where  Z  *  C  .  The  variables  are  taken  to  be  X^  and  Z  in  Eq.  (21),  to  avoid 
a  singularity  at  Z^.  At  Z^,  Eq.  (21a)  becomes 
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G  F(0,X  ) 
c  *  e 


(1  -  2Z2  -  2Zl  G.) 


(22) 


which  is  independent  of  R.  A  numerical  integration  of  Eq.  (21)  is  generally 
required. 

The  intensity  distribution  is  obtained  from 


tt  Av  I  +  1+R  .  /  .  \ 

TTS-g-r - - - -  (w  ‘  Z)  +  RH(Xf,Xa) 

F(X’V  G  F(0,X  )  \2Z  /  f  6 

c  e 


(23) 


The  output  power  to  station  Z  equals,  from  Eq.  (14), 


2?d/2  P 


^Z  -  f  Z 3  -  Gc  F(0,Xe)  H(Xf,Xe)J 


-  2  Gc  F(0,Xe)  Jz  Z  H(XffXe)  dZ 


(24) 


Numerical  results,  presented  in  Tables  I  through  III  and  Figs.  5  through  7, 
are  discussed  in  Section  III. 

3.  LIMIT  R  »  1,  RG  >1 
’  c 

The  limit  R  >  1,  R  G£  >  1  is  now  considered.  It  is  clear  from  physical 

considerations  that  Xe  decreases  as  R  and  G  increase.  The  present  limits 

2  c 

ensure  that  terms  of  order  X^  can  be  neglected  compared  to  1. 
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Table  I.  Station  at  which  lasing  is  initiated 


Table  III.  Dovmseream  station  at  which  Xj  -  0  and 
corresponding  output  power 


Xe  - 

0.001 

*e 

-  1.0 

X« 

-  2.0 

I  J/2-  (.) 

P 

r 

PL 

7. 

PL 

Pt 

2  T)  c 

Pe.s 

Pe,s 

zl 

Pe.s 

3.53G"3 

1 

0.978 

0.727 

0.977 

0.722 

0.975 

0.713 

10 

0.786 

0.964 

0.786 

0.957 

0.783 

0.949 

100 

0.724 

0.985 

0.724 

0.981 

0.720 

0.971 

1000 

0.710 

0.988 

0.710 

0.984 

0.707 

0.974 

m 

0.705 

0.989 

0.705 

0.985 

0.703 

0.975 

7.060-3 

i 

0.977 

0.706 

0.977 

0.696 

0.972 

0.677 

10 

0.789 

0.946 

0.790 

0.936 

0.785 

0.917 

100 

0.724 

0.973 

0.725 

0.965 

0.718 

0.945 

1000 

0.710 

0.976 

0.709 

0.969 

0.704 

0.948 

0.704 

0.978 

0.701 

0.970 

0.699 

0.949 

1.765*2 

l 

0.974 

0.650 

0.971 

0.628 

0.960 

0.585 

10 

0.795 

0.898 

0.795 

0.875 

0.784 

0.827 

100 

0.725 

0.936 

0.723 

0.916 

0.711 

0.867 

1000 

0.706 

0.942 

0.704 

0.923 

0.694 

0.872 

- 

0.698 

0.944 

0.695 

0.925 

0.687 

0.874 

3.530"2 

1 

0.966 

0.572 

0.960 

0.531 

0.937 

0.454 

10 

0.798 

0.825 

0.796 

0.783 

0.777 

0.689 

100 

0.720 

0.877 

0.716 

0.839 

0.696 

0.741 

1000 

0.700 

0.886 

0.693 

0.848 

0.677 

0.748 

m 

0.690 

0.889 

0.684 

0.851 

0.668 

0.750 

7.060”2 

l 

0.949 

0.441 

0.934 

0.374 

0.882 

0.244 

10 

0.796 

0.693 

0.790 

0.615 

0.753 

0.437 

100 

0.710 

0.764 

0.701 

0.690 

0.665 

0.501 

1000 

0.684 

0.776 

0.674 

0.703 

0.641 

0.509 

• 

0.673 

0.779 

0.661 

0.706 

0.632 

0.512 

‘Superscript  denotes 

exponent 

of  ten 

(••8*  t  3« 

530“3  -  3. 

530  x  10 

-3>. 
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FREQUENCY  WIDTH,  X, 


Figure  6(a).  Numerical  solution  of  Eqs.  (10),  (21),  and  (24)  for 
caae^of  laminar  flame  sheet  R  »  1,  X_  •  0.001,  and 
10  •  0.353,  1.765,  7.060;  frequency  width,  Xj 
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Numerical  solution  of  Eqs.  (10),  (21),  and  (24)  for 
case^of  laminar  flame  sheet  R  •  1,  X£  *  0.001,  and 
10Z  G  ■  0.353,  1.765,  7.060;  net  inversion,  AN/AN. 


OUTPUT  POWER,  P/P 


i 


Figure  6(c).  Numerical  solution  of  Eqs.  (10),  (21),  and  (24)  for 
caae^of  laminar  flame  sheet  R  ■  1,  Xg  «  0.001,  and 
10Z  G£  ■  0.353,  1.765,  7.060;  output  power,  pe/pe>s 
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j  ?5  i  - EXACT  SOLUTION  [EQ.I21I] 

- LIMIT  fl»1,GrR  =  0111  [EQ.I29I] 
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STREAMWISE  DISTANCE,  Z 

Figure  7(a).  Numerical  solutions  of  Eqs.  (10),  (21),  and  (24)  for 
.  caae^of  laminar  flame  sheet  R  “  100,  Xfi  =  0.001,  and 
10Z  G  *  0.353,  1.765,  7.060,  including  comparison 


SOLUTION  [EQS.  1211, 1101) 


l-  NV/NV  'NOISdTANI  13N 


rical  solution  of  Eqs.  (10),  (21),  and  (24)  for 
of  laminar  flame  sheet  R  *  100,  Xe  =  0.001,  and 
G  -  0.353,  1.765,  7.060,  including  comparison 
analytic  solution  in  limit  R  »  I,  R  C  ■  0(1); 
inversion,  AN/AN 


Substitution  of  the  identity 


1  -  H(Xf,Xe)  F(Xf,Xe) 


-X2 

^T7I  xf  [*  + 


(25) 


3 

into  Eq.  (21a)  indicates  RX.  is  of  the  order  1  in  the  present  approximation. 

£  2 

Neglect  of  terms  of  the  order  X^  compared  with  1  in  Eq.  (21a)  Indicates  that 
the  derivative  terra  is  negligible  therein.  Equation  (21a)  becomes 


t  ,  1/2  YT 

X3  -  -  *  e  6 
f  4  ~ 

G  R 
c 


Other  quantities  of  interest  are 


-X 

6  2 

AN  '  6  Xf 

ANX  *  F(0,Xe) 


3ir 


172  e  X*  X3  +  0  (xj) 


IT  Av,  G  I 
h  c 


-X 

G  R  e  6  . 

f(o,x  )  ^Xf 
6 


X2)  [1  +  0  (x2)] 


2  4/2  %  * 1  •  “2  • 2  5C  z[i  +  0(xt)] 


2  «D/2  p  -  jz  -  f  Z3  -  Se  z2  [1+°(Xf)]j 


(26a) 


(26b) 


(26c) 


(26d) 


(26e) 


40 


The  station  Zft  at  which  dP/dZ  ■  0  is 


Ze  -  [(2  +  G^)1/2  -  Gj/2 


(27) 


Net  output  power  Is  obtained  by  using  Zg  as  the  upper  limit  in  Eq.  (26e). 

Results  for  a  saturated  laser  can  be  obtained  by  considering  the  limit 

R  +  •,  G  ♦  0,  RG  »  1.  It  is  found  that 
c  c 


Z 


i 


0 


(28a) 


P 

e,s 


(28b) 


where  P„  denotes  net  saturated  output  power  and  is  a  convenient  reference. 
e»S 

In  the  present  limit,  the  dependent  variables  X-,  AN/AN.,  and  G  I  depend 

~  1/2  r  1  c~ 
on  Z,  Xg,  and  GcR.  The  parameters  Z±t  Z&,  and  ?D  Pg  depend  only  on  Gc«  The 

medium  acts  like  a  homogeneously  broadened  medium,  i.e.,  X^  *  0,  and  lasing 

occurs  at  line  center.  Equations  (26b)  and  (26c)  are  the  same  as  those 

obtained  for  a  homogeneous  medium  with  lasing  at  line  center  and  with  zero 

power  line  center  gain  reduced  by  the  factor  F(0,Xe)  to  account  for  the  effect 

of  spreading. 


4.  LIMIT  R  »  1,  RGc  -  0(1) 

The  present  limit  R  »  1,  RGc  ■.  (1)  implies  Gc  «  1  and,  therefore, 

implies  a  saturated  laser.  Moreover,  because  of  the  conflicting  influence  of 

R  and  G  on  X.,  the  limit  RG  ■  0(1)  ensures  that  Xr  is  of  order  1. 
c  r  *  c  i 

If  terms  of  the  order  Gc  and  R“*  are  neglected,  the  derivative  term  in 

eq.  (21a)  is  again  negligible,  and  Eq.  (21a)  becomes 
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(29a) 


Z  -  ((2  +  $2)1/2-  SJ/2 


where 


0 


Gc  R  [1  -  F(Xf,Xe)  H(Xf,Xe)]  F(0,Xe)/F(Xf,Xe) 


(29b) 


An  implicit  solution  for  X^  in  terms  of  Z  is  provided  in  Eq.  (29).  It  is 
found  that  X^  decreases  monatomically  as  Z  increases  from  Z  ■  0  to  Z  -1//2. 
In  the  vicinity  of  Z  ■  0  and  Z  ■  1//2,  respectively, 


Xj  -  (-1)  In  [2F(0,X  )  G  RZ] 
t  e  c 


(30a) 


3ir 


1/2 


ti  (?k  -  j  t1 + 0(x2f)] 


(30b) 


Thus,  Xj  ♦  •  as  Z  +  0,  and  Xj  +  0  as  Z  ♦  1//T.  The  derivative  dXf/dZ  becomes 
infinite  as  Z  +  0  and  as  Z  ♦  1//2,  and,  therefore,  should  be  retained  in  Eq. 
(21a)  when  these  regions  are  considered.  Lasing  is  initiated  at 
station  Z^  *  Gc  «  1.  The  boundary  condition  Xf  *  0  at  is  not  satisfied 
because  of  neglect  of  the  derivative  term  in  Eq.  (21).  Lasing  intensity  is 
found  from 


v  Av  G  I  -  G 
he  c 


R 


'  F(X,Xe) 
F(XflXe) 


(31) 
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The  Intensity  is  proportional  to  zero  power  lineshape,  and  T!  ♦  0  as  X  +  Xj, 
Since  <?c  «  1,  the  laser  is  saturated  and  output  power  is  obtained  from 
Eq*.  (28a)  and  (28b).  Equations  (28),  (30),  and  (31)  indicate  that  X^,  dP/dZ, 
sad  I  all  become  zero  at  Z  •  1//2.  In  more  general  cases,  neither  Xf  nor  I 
is  zero  at  the  station  where  dP/dZ  ■  0  (i.e.,  at  Ze). 

In  the  present  limit,  the  variables  Xf,  AN/dN,,  and  G  I  depend  on  Z,  X_, 

-  1/2  1  C 
and  CcR.  The  quantities  Z^,  Zfi,  and  P  ,  however,  are  constants  in  the 

limit  G  «  1. 
c 

An  improved  estimate  for  output  power  can  be  obtained  by  using  a  mean 
value  for  H(X£,Xg)  in  the  integral  of  Eq.  (24).  The  result  can  be  expressed 


where  the  error  term  tends  to  become  small  as  Z  ♦  Zg  (e.g.,  Fig.  7).  Net  out¬ 
put  power  is  obtained  by  utilizing  Ze  from  Eq.  (27)  as  the  upper  limit  in 
Eq.  (32).  Equation  (32)  is  equivalent  to  Eq.  (26e),  and  both  provide  accurate 
estimates  of  laser  output  power  for  R  »  1. 


III.  RESULTS  AND  DISCUSSION 

In  a  typical  cv  chemical  laser,"  R  ■  0(10)  -  0(100)  and  G  -  0(0.01). 

c 

Hence,  the  Halting  solutions  In  Sections  II. C. 3  and  II. C. 4  Include  aost  cases 
of  practical  Interest.  The  intensity  I  is  positive  for  all  values  of  C  and 
Xf,  and,  therefore,  these  Halting  solutions  are  physically  realistic.  The 
negative  values  of  I  encountered  at  lower  values  of  R  can  be  avoided  if  the 
lineshape  in  region  |X/Xfl  >  1  is  not  specified  a  priori.  The  regime  where 
negative  values  of  I  occur  Is  of  less  interest,  and  realistic  solutions  for 
this  regime  are  not  pursued  herein. 

Numerical  results  have  been  obtained  for  the  case  of  a  laminar 

flame  sheet  [Eqs.  (10),  (21),  and  (24)].  Parameters  in  the  range  0.00353  < 

1  /2 

( 2/ir  )Cjj  G£  <  0.0706,  0.001  <  Xg  <  2.5,  and  1  <  R  <  •  were  considered.  The 

first  of  these  inequalities  corresponds  to  0.01  <  (g  )  /g  <0.2 

c  av  mzp 

[Eq.  (19)].  Results  are  presented  in  Tables  I  through  III  and  Figs.  5 
through  7. 

In  Table  I,  the  values  of  Z^  at  which  lasing  is  Initiated  are  Indicated. 

These  are  Independent  of  R  and  increase  as  X@  and  Gc  Increase.  In  Table  II  is 

presented  the  station  Ze  at  which  dP/dZ  ■  0  and  the  corresponding  output  power 

P.  normalized  to  the  saturated  laser  output  power  P_  _.  The  ratio  P./P.  _ 

“  6)8  “  “  )8 

decreases  as  X£  and  Gc  increase  and  is  relatively  insensitive  to  R  for  R  >  0 
(10);  therefore,  Eq.  (26e)  provides  the  output  power  for  cases  of  practical 
interest.'  The  results  for  R  ■  100  are  also  indicated  in  Fig.  5.  The  down¬ 
stream  station  at  which  X^  *  0  is  denoted  Z^.  Values  of  Z^  and  the  corre¬ 
sponding  output  power  P^  are  indicated  in  Table  III.  The  values  of  P^  depart 
significantly  from  the  corresponding  values  of  P£  in  Table  II  only  for  R  -  1. 
The  difference  is  the  result  of  negative  values  of  I.  Thus,  the  data  of 
Table  III  Illustrate  the  nature  of  the  present  solution  for  R  ■  0(1)  and  are 
primarily  of  academic  Interest. 
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The  variation  of  Xf,  AN,  and  P  with  Z  Is  Indicated  in  Fig.  6  for  R  *  1.0, 
Xe  ■  0.001,  and  several  values  of  G£.  Results  from  the  numerical  integration 
of  Eq.  (21)  for  R  ”  100  are  compared  with  results  from  the  limiting  solution 
of  Section  II. C. 4  in  Fig.  7.  The  agreement  is  good  except  for  points  in  the 
vicinity  of  and  Ze  where  poor  agreement  is  expected  because  of  neglect  of 
the  derivative  term  in  the  limiting  solution. 
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IV.  CONCLUDING  REMARKS 


In  Table  II  laser  output  power  is  demonstrated  to  depend  primarily 

on  G  and  and  to  be  relatively  independent  of  variations  in  R  for  the 

practical  range  R  >  0(10).  Hence,  solutions  in  the  limit  R  ♦  •  provide  useful 

estimates  for  the  effect  of  gain  saturation  G  and  source  flow  on  the  out- 

c  .  e 

put  power  from  multiple  longitudinal  mode  cw  chemical  lasers. 

A  number  of  complex  mult itrans it ion-multichemical  reaction  numerical 
codes  have  been  developed  (e.g.,  Ref.  8)  that  can  be  used  to  evaluate  cw 
chemical  laser  performance  with  the  assumption  of  a  single  longitudinal  mode 
at  line  center  for  each  lasing  transition.  In  the  present  model,  X^  ♦  0  as 
R  >  •,  i.e. ,  the  limit  R  <*  corresponds  to  a  single  longitudinal  mode  at  line 
center.  Therefore,  numerical  codes  in  which  this  assumption  is  made  provide 
reasonable  estimates  for  the  effect  of  gain  saturation  and  source  flow  on 
multiple  longitudinal  mode  cw  chemical  laser  output  power  if  the  correct  zero 
power  line  center  gain  is  used  [Eq.  (5)],  and  the  conditions  in  Eqs.  (la)  and 
(lb)  are  satisfied. 

It  is  also  seen,  from  Fig.  7,  that  the  limit  solution  R  »  1  and  G  R  - 

c 

0(1)  provides  reasonably  accurate  simple-closed-form  analytic  expressions  for 
laser  performance  in  the  regime  of  interest. 
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SYMBOLS 


D(  ) 
F(X,Xe) 
G(v),  Gc 
G 

c 

g(v) 

8c»®mzp 

H(Xf,Xe) 

Ij.ICv) 

J 

kcd*kcr 

L 

N.Nj.Nj.N,. 

H(v),Nv 

an,anv 

n,n(v) 


R 


Most  probable  particle  thermal  speed 
Dawson  integral,  Eq.  (17) 

Zero  power  lineshape,  Eq.  (5) 

Normalized  gain  and  threshold  gain,  Eqs.  (4a)  and  (7a) 
Threshold  gain  parameter,  Eq.  (20) 

Small  signal  gain,  Eq.  (4) 

Threshold  gain,  maximum  zero  power  gain,  Eqs.  (7b)  and  (19) 

Net  inversion  AN/AN  ,  Eq.  (10) 

Vi 

Lasing  intensity,  Eq.  (11) 

Longitudinal  mode  index 

Collisional  deactivation  and  cross  relaxation  (molecular 
collision)  rates,  sec-1 

Mirror  separation 

Normalized  population  density,  Eqs.  (2),  (11),  and  (16) 

Normalized  population  in  frequency  interval  v  to  v  +  dv, 

Eq.  (2) 

Normalized  population  difference,  Eq.  (2) 

Population  density,  Eq.  (2) 

Number  of  semichannels,  Fig.  1 

Output  power  liberated  up  to  station  c,  ;e»  Eq.  (11) 

Output  power  for  saturated  laser,  Eq.  (28) 

Reciprocal  of  Doppler  width,  Eq.  (2) 

kcr^kcd 


51 


r 


w 

X 

Xe 

*f 


x,y 


*0 

yf  (x) 


c,cD 


V 

Av  ,Av 
c 

a 


d-Avh 


Velocity  In  x  direction.  Fig.  1 

Transverse  velocity  at  edge  of  source  flow.  Fig.  3 
Channel  semiwidth.  Fig.  1 
Normalized  frequency,  Eq.  (5) 

Source  flow  parameter,  Vg/ a 

Value  of  X  corresponding  to  largest  lasing  frequency 
Streamwise  and  lateral  distance.  Fig.  1 
Characteristic  diffusion  distance.  Fig.  1 
Flame  sheet  location,  Fig.  1 

Normalized  streamwise  distance,  kcdx/u,  kC(jxD/u 
Frequency,  sec-* 

Characteristic  frequencies,  Eq.  (1) 

Stimulated  emission  cross  section,  Eq.  (4) 


Subscripts 


o 

1,2 

av 

c 

e 

f 

i 

i 


Line  center  value 
Levels  1  and  2 
Average  value,  Eq.  (4) 

Related  to  cavity 

End  of  positive  lasing  region  (station  at  which  dP/dc  -  0) 
Associated  with  final  (highest)  lasing  frequency 
Value  at  start  of  lasing 
Longitudinal  mode  index 
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o 


Downstream  station  at  which  Xj  - 
Pertaining  to  frequency  v 
Pertaining  to  highest  lasing  frequency 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  la  conducting  exp ar¬ 
isen  tal  and  theoretical  Investigations  necessary  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  military  space  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  in  the  nation's  rapidly  developing 
space  systems.  Expertise  In  the  latest  scientific  developments  Is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophyslcs  Laboratory :  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high-temperature  thermomechanics,  gas  kinetics  and  radiation;  research 
In  environmental  chemistry  and  contamination;  cw  and  pulsed  chemical  laser 
development  Including  chemical  kinetics,  spectroscopy,  optical  reeonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuus  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory :  Microelectronics,  CaAs  lov-noiae  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  mllllmeter-vave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation  performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  Intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  In 
space  environment;  materials  performance  In  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  In  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  Ionospheric  physics,  radiation 
from  the  etmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  airglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  In  the  magnetosphere;  solar  physics.  Infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere.  Ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  In  space  on  space  systems. 


